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Abstract 

We study theoretically two effects related to impurity depinning that are common for ID Luttinger liquid (LL) and 
linear-chain CDW conductors. First, we consider the electron transport through a single impurity in a ID conducting 
channel and study a new regime of conduction related to LL sliding at voltage above a threshold one. The DC current 
in this regime is accompanied by oscillations with frequency / = I/e. This resembles the CDW depinning in linear- 
chain conductors, the Josephson effect, and the Coulomb blockade. Second, we found that strong pinning of the LL 
by impurities leads to a magnetic field dependence of the low-temperature specific heat similar to that observed 
experimentally in CDW compounds. We interpret this in favor of possibility of formation of the LL in linear-chain 
compounds. 
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1. Introduction 

In contrast to 3D electronic systems where elec- 
trons form, typically, a Fermi liquid, electrons in 
ID systems are known to form a Luttinger liq- 
uid (LL) [T]. In case of the long-range Coulomb 
interaction the LL can be considered as a ID ver- 
sion of the Wigner Crystal [2]. So the LL can be 
viewed as a simple prototype of an electronic crys- 
tal. A common feature of the LL with repulsive 
inter-electronic interaction and of other electronic 
crystals, like Wigner Crystal or charge-density wave 
(CDW), is their pinning by impurities. It is well- 
known that impurities in ID systems with repulsive 
inter-electronic interaction form effectively large 
barriers and strongly suppress the current that 
leads to a power-law dependence of conductivity 
on voltage and/or on temperature [3]. This effect 
can be described in terms of macroscopic tunnelling 
between different minima of a washboard poten- 



tial describing interaction with an impurity [T]. In 
Sec. [2] we study a new regime of conduction in ID 
LL that corresponds to depinning of LL and its slid- 
ing at higher voltages corresponding to a large tilt 
of the washboard potential permitting depinning 
of the LL from the impurity. The DC current / in 
this regime is accompanied by coherent oscillations 
with frequency / = I/e. The effect is similar to the 
well-known electric field depinning accompanied by 
generation of the narrow-band noise in linear-chain 
CDW compounds [1] . 

The concept of the Luttinger liquid can be ap- 
plied to quasi- ID CDW compounds as well. It 
was shown theoretically |5|l6j that formation of the 
CDW in quasi- ID materials with conducting chains 
coupled by elastic forces can be described in the 
frame of the LL if the interaction between the elec- 
trons is strong enough. And, in fact, the electronic 
interaction in these materials is, typically quite 
strong. There are also experimental indications to 
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LL-like behavior of linear-chain compounds: a tran- 
sition to a low temperature state characterized by 
LL-like conductivity was detected in focused-ion 
beam processed crystals of TaSa and NbSea [7] . In 
order to account for such behavior the possibility of 
stabilization of the LL state by defects in quasi-lD 
metals was suggested [5]. In general, whereas the 
main properties of quasi- ID conductors at higher 
temperatures are well understood, at low tempera- 
tures they demonstrate many intriguing properties 
which are not yet explained convincingly, e.g., un- 
expected strong magnetic field dependence of the 
specific heat in non-magnetic CDW materials TaSa 
and blue bronze 9 lO'll]. In Sec. [3] we study the 
specific heat of pinned CDW described in terms of 
the LL, and find a magnetic field dependence of the 
heat capacity similar to that observed in the CDW 
materials. This dependence originates from strong 
pinning of the LL by impurities, and on partial 
violation of spin-charge separation by impurities 
resulting in possibility of localization of an extra 
electron spin in segments between impurities. 

2. Depinning of the Luttinger liquid from an 
impurity in ID conducting channel 

The LL Hamiltonian describing interaction with 
impurities at a; = a;^ reads [I] 

Himp = -^Yl J dxW cos{V2^p{t,x) ~ 2kFXi) 

i 

X C0sV2^a{t,x) 5{x - Xi), (1) 

where with v = p,a are the bosonic displacement 
fields in charge (spin) sectors, and W is related to 
2A;p'-component of the impurity potential. 

In this section we study the response of the repul- 
sive LL with single impurity at a; = to DC volt- 
age. For brevity we consider first the spinless LL us- 
ing the Tomonaga-Luttinger (TL) model with short 
range interaction between electrons. Such interac- 
tion describes gated quantum wires where the long- 
range part of the interaction is screened by 3D gate 
electrodes. At the end we discuss essential modifi- 
cations induced by the long-range Coulomb interac- 
tion and spin. 

Periodicity of the interaction potential is associ- 
ated with Friedel oscillations induced by impurity. 
Note that impurity pinning of the CDW is also re- 
lated to Friedel oscillations [12]. The particle den- 



cos {2kpx — 2$p), where the first term describes 
smooth variations of the particle density. Fluctua- 
tions of $ make the expectation value of the rapidly 
oscillating second term in a free LL equal to zero. 
However, near the impurity the second term survives 
resulting in Friedel oscillations, i. e., in 2kp modu- 
lation of charge density around an impurity. As long 
as the current operator in the LL is given by relation 
/ = {e/7T)dt^p the current flow through the impu- 
rity means an increase of $p with time, which cor- 
responds to a shift of the Friedel oscillations. So the 
Friedel oscillations is an obstacle for the current flow 
through an impurity. The well-known power-law I-V 
curves |ll3j are induced by tunneling between min- 
ima of washboard potential ([1]) slightly inclined by 
an external voltage. This resembles Josephson junc- 
tions where fluctuations result in a flnite voltage 
drop at currents below the Josephson critical cur- 
rent, I < Ic (for a review see Ref. However, in 
superconducting junctions at I > Ic the Josephson 
oscillations start, and this corresponds to an increase 
of the superconducting phase difference with time in 
the washboard potential inclined to an amount ex- 
ceeding the critical value. Below we show that a sim- 
ilar regime must occur in the LL with an impurity 
when the applied voltage exceeds a threshold slope 
of the washboard potential at which the system can 
roll out from the minimum. Above the threshold the 
current is larger than the tunneling current in the 
sub-threshold regime corresponding to power-law I- 
V curves and is accompanied by oscillations with the 
washboard frequency / = I /e. 

The spinful Hamiltonian includes the impurity 
term ([T]), the standard TL Hamiltonian of interact- 
ing electrons 



H' 



TL 



E 



dx 



hlTVF 



,(2) 



and the term with an external electric field E 



where Kp is the LL parameter measuring the 
strength of the interaction: Kp < 1 for repulsive 
and Kp > 1 for attractive interaction. For a spin- 
independent electronic repulsion K^ = 1. 

Now we return to a spinless case and till the end 
of this section drop the subscript at $p. Commut- 
ing field operators with the Hamiltonian we derive 
equation of motion for the Heisenberg operator $ 



sity operator in LL is expressed as p 



1 p 

TV dx 



{vldl 



(3) 
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2W sm2^o{t)S{x) - E 



where $o(i) = ^{t^x — 0), Vp — vp/Kp is the 
plasmon velocity, and 7 describes small damping 
induced by coupling of electrons to a dissipative 
bosonic bath (to phonons or to density fluctua- 
tions in a metallic gate |14I15| V Using Eq. ([3]) with 
boundary conditions at the contacts one can ex- 
press $(t,a;) in terms of its value at the impurity 
site. We apply boundary conditions _16J derived for 
a wire adiabatically connected to ideal Fermi-liquid 
reservoirs at x = ±L/2 with voltage difference V . 
We consider a long wire, L ^ w_f/7, and neglect 
reflection of current pulses generated by the impu- 
rity from the contacts. Then we obtain equation of 
motion for $o(t) at the impurity position supple- 
mented by relation between time averaged values 

00 

5t$o(t) + J dheZ{t-h)sm2^o{h) = ^l (4) 


/ = GoiV Zit) = W,Kp [ f^e--* 

J 211 V ^ + *7 

where / = e^i?Wi?/7r?i7 is the DC current, Go = e^/h 
is the conductance quantum per spin orientation, 
and 



= 2W^(sin2$o(t))t 



(5) 



is the DC voltage drop at the impurity. 

As $0(0 is operator, it is not easy to solve 
Eq. (U]). However, in the limit of strong inter-electron 
interaction {Kp 0) when fluctuations of $o(i) at 
the impurity are small, '^^{t) can be treated as an 
ordinary function. So we consider first the simple 
case Kp — > when fluctuations can be neglected 
and Eq. (|4|) can be easily solved analytically. When 
V < Vt = 2W the solution is stationary, 2$o — 
arcsin (C//Vt), with zero current / — 0. At V > 
Vt the solution is oscillatory with fundamental fre- 
quency f — I/e 



TTp/e 



+ {GoVtKp)^ + GoVtKp sin 



(6) 



Eq. ^ determines the current at the impurity site, 
I{t,x = 0) = e<i'o(i)/7r. Current at the clean part of 
the channel is equal to I{t, x) — edt^o{t— \x\/vp) at 
|a;| ^ w_f/7, and /(i, a;) = / at large distances from 
the impurity. 

For the DC voltage drop at the impurity we obtain 



P + {GoVrKpy ~ I 



/(GoKp). 



(7) 



Thus the oscillatory regime starts a.t V > Vt and 
the amplitude of the oscillations decreases as voltage 
increases decaying at large voltages V ^ Vt- 

Now we consider the case of finite values of 
Kp. Mean square fiuctuations ((5$^) can be cal- 
culated from Keldysh Green's function = 
-i{{d^(t),6^{t')}+) at t = t'. This function can 
be expressed via retarded and advanced Green's 
functions 

D^'^^\t,t') ^ ±ie[±{t - t')]([s^{t),s^{t')]^) 

by relation i:>^'(i,i') = D^{t,t')f{t')-f{t)D'^{t,t'), 
where Fourier transform of / is related to the dis- 
tribution function of bosonic excitations N(uj), 
f{u!) = 1 + 2N{uj). In the equilibrium state N{uj) is 
the Planck distribution function. At low tempera- 
tures, smaller than all characteristic energies of the 
system, one can neglect contribution of thermally 
excited excitations, then f{io) = sign((jj) and 
can be expressed via the retarded and advanced 
functions. Note that acting in this way we neglect 
the effect of non-equilibrium distribution of bosonic 
excitations. 

Now we derive equations of motion for the re- 
tarded and advanced Greens functions of fluctua- 
tions. This can be done in a standard way multi- 
plying Eq. [3]by ^{t',x') from the left and from the 
right, and combining them in order to obtain cor- 
responding Green's function after averaging. Then 
using the Fourier transformation we express Greens 
functions at the impurity site and get rid of the co- 
ordinate dependence of Green's function. After that 
using the statistically averaged equation (jH) we sub- 
tract expectation values and use the self-consistent 
harmonic approximation 

sin2(5<l -> (cos2<5$)254 = e-^^^^''h5^, 

and arrive, flnally, to close equations of motion for 
jjR{A) ^ As we need this equations at frequencies 
larger than the small damping constant we neglect 
damping 7 and employ Z{t) — WKpS{t — ti). Then 
the equation for reads 



2e 



ttK 



dtD^ + -WKpC{t)D" 2 
G{t) = cos2$(ti)(cos2(5<l(ti)). 



(8) 



This equation and the similar equation for the ad- 
vanced Green's function can be easily solved analyt- 
ically. The solutions are 
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This gives us for {5^"^) = ^D^{t, t) 

OC 

2 J tl 



Eq. pO)) must be solved self-consistently with the 
statistically averaged equation (|4]) in order to find 
((5$^) as function of cos2<i>o- First, we calculate the 
threshold voltage. In the sub-threshold region C 
does not depend on time. Then calculating integrals 
in pn)) and using definition of C ([5]) we obtain 

/5$2\ ^ ^ 1 ^ 

^ ' 2{l -Kp) 2KpeWcos2^n 

where K ^ ppv is a, large cut-off energy. In accor- 
dance with our previous statement we found that 
fluctuations vanish at Kp 0. Substituting (fTTjl in 
([5]) we obtain equation for $o in the sub-threshold 
regime 

^^^/2KpeWcos2<^o\^-'''' . 

i — A sin2$o = ^^. (12) 

It has solution a,t V < Vr with 

v,, = w(B£l^)'^yT^^ (13) 

We see from ([T3l) that the finite threshold voltage ex- 
ists provided Kp < 1, and this is in accordance with 
the condition that the impurity is a relevant per- 
turbation in case of repulsive inter-electronic inter- 
action. Note, however, that fluctuations were taken 
into account in the self-consistent harmonic approx- 
imation that does not work well at Kp — )■ 1. 

It is not simple to take into account fluctuations 
analytically in a general time-dependent case, but 
this can be done easily near the threshold value when 
the DC current is small, / ^ GoVrKp. In this case 
we find ^ and ([7]) again, but with different value 
of the threshold voltage Vr given by (fT3| . 

So far we considered the spinless LL. In the spin- 
ful LL the spin-charge separation is violated at the 
impurity. This leads to modification of the results, 
in particular, of the threshold voltage. For a spin- 
independent electronic repulsion Eq. (|13p is substi- 
tuted by 

v, = „,(H£E)^^-7^^/^^. 



Note also that electron transport current is accom- 
panied by pulses of spin current generated at the 
impurity. 

We considered also the case of the long-range 
Coulomb repulsion, and find that it does not re- 
sult in important qualitative difference. One of the 
distinctions is the smaller role of fluctuations. An 
estimate for the threshold field in this case gives 
with the logarithmic accuracy 

where d is the wire diameter, — -4^ — is the 
dimensionless parameter measuring the strength of 
the Coulomb interaction, and e is a background di- 
electric constant. 

In conclusion, we found that above the threshold 
voltage Vt the current through an impurity gener- 
ates coherent oscillations with the fundamental fre- 
quency / — I /e. If in addition there is also an AC 
applied voltage with frequency /o then an analog of 
the resonance Shapiro steps observed in Josephson 
junctions will appear on the I-V curves. In our prob- 
lem these are the steps of a constant voltage, the 
fundamental step being located at the current value 
/ — e/o. Characteristic frequencies of the oscilla- 
tions are determined by the strength of the impurity 
potential. In semiconducting quantum wires typical 
values of the impurity potential can be of the order 
of several meV, so depending on the strength of the 
electronic repulsion the frequency may fall into gi- 
gahertz or terahertz frequency region. Direct appli- 
cation of our results to real systems are limited by 
voltages smaller than distances to other electronic 
subbands. The results can be modified also due to 
different coupling of the ID system to 3D environ- 
ment. 

3. Magnetic field dependence of 
lov^r-temperature specific heat of quasi-lD 
conductors 

Now we take into account formation of the CDW 
due to electron-phonon coupling. We describe free 
phonons and electron-phonon coupling are de- 
scribed in a standard way, free phonons being pre- 
sented by an elastic displacement field ip with a 
harmonic action. We need to keep only slowly vary- 
ing part of the 2fci? component of ip leading to the 
CDW transition, therefore, in most of the cases the 
weak dependence of the phonon frequency Ljph{c\) 



4 



on longitudinal component of the momentum, , 
can be neglected. Then we bozonize the action 
(confer, e.g., [17]) and integrate out phonon degrees 
of freedom. This leaves the effective inter-electron 
interaction term 

•^int = (~) X/ / dxdTdTiD^{T - Ti,n- m) 

X cos[A/2i>p(T, X, n) - %/2$p(ti, a:, ni)] 

X cos v^^ctIt, n) cos \/2$o-(ti, a;. Til) (14) 

where g is the electron-phonon coupling constant, 
a is a small cut-off distance, and is the free 
phonon propagator. This equation differs from that 
of Ref. [17j by summation over the chain numbers n. 

With such action formation of the CDW was stud- 
ied in [6]. In ou problem we must add the impuri- 
ties ([T]). We concentrate on a simple case of very 
strong inter-electronic interaction {Kp <C 1) result- 
ing in strong pinning and small fluctuations of the 
fields <I>^ around classical trajectories corresponding 
to the minima of the total action. Fot the case of 
the strong pinning we must find classical trajecto- 
ries satisfying the boundary conditions at the impu- 
rity positions. Boundary conditions that deliver the 
minimum of the impurity action are 

p{t,Xi) — 2kFXi = nn, V2^cr{t, Xi) = tott 

where (n + m) must be an even integer. The let- 
ter condition means that the number of electrons in 
segments between the impurities is an integer. 

Consider first, for illustration, the case without 
the CDW when the electron-phonon coupling is ab- 
sent. In this case the system is effectively broken into 
independent segments with bounded LL. The energy 
spectrum of the segments consists of zero modes and 
excitations, the energy of the zero modes being de- 
pendent on number of extra electrons at a segment 
confined by impurities (confer [8]). The lowest en- 
ergy of the zero-mode state corresponds to number 
of extra electrons 5n = Q ov 5n — ±1 depending on 
the values of the phases at the impurity. The states 
with one extra or one missing electron behave sim- 
ilar to magnetic impurity since they can have spins 
of both directions, and this leads to magnetic field 
dependence of the specific heat. Exited states are 
separated from the zero mode states by energies of 
the order of v^/lKp, where Vi, = vp/K^, v ~ p, ct, 
and I is a typical distance between the impurities, so 
the excited states can be ignored if the temperature 
is low enough. 



A similar picture takes place in the presence of the 
CDW, the main difference being related to a coor- 
dinate dependence of and $cr between the impu- 
rities since the phase perturbations will be localized 
near the impurities. An exact solution for the zero 
mode phase distribution depends on details of the 
crystalline structure and phonon dispersion. There- 
fore, we concentrate on the qualitative features and 
consider a simple model resembling approaches used 
to describe interaction of the CDW with impurity in 
Refs. |18|19j (the spin-density wave case was treated 
similarly in ppj). In these papers phase variation 
around an impurity was described by static sine- 
Gordon equation. Our case is somewhat different 
because we need to solve coupled equations for two 
phases $p and Note that the details of phase 
variations are not important for our problem, fur- 
thermore, we will find that the magnetic field de- 
pendence of the specific heat is rather universal and 
does not depend on these details and on energies of 
the zero-mode states. So we use an oversimplified 
model analogous to that of Refs. |18|19j . Namely, 
having in mind small impurity density we consider 
isolated impurities and neglect in action (|14p vari- 
ations of the phases at the chains that does not 
contain the impurity. Matching the zero mode so- 
lutions of the equations for the phase variations at 
the impurity we found that the smallest energy can 
be achieved by the zero-mode states with number 
of extra electrons localized near the impurity 8n — 
0, ±1. These states correspond to variations of the 
phases 5$p and 5$cr around the impurity equal to 
(0,0) and (±7r/\/2, ±7r/\/2), respectively. The for- 
mer state has energy iJp.o = 7^(1 — cos^i) in the 
charge sector and Eafs = in the spin sector, where 
M is an energy of the order of the CDW gap, and 
— kpXi. The latter states describing charge ±e 
and spin ± 1 localized near the impurity have energy 
Ep,i — -^{2 — cos0i — sin0j) in the charge sector 
and Ea.i ^ M <^ Ep,i in the spin sector. In the 
presence of magnetic field the energy of the states 
with ±1 extra spin must be shifted by ±iibH, where 
/is is of the order of the Bohr magneton of an elec- 
tron. The lowest excited states in the charge sector 
are separated by the energy gaps of the order of a 
phonon frequency ujp ^ s±/d where is the sound 
velocity in the direction perpendicular to the chains, 
and d is an inter-chain distance. In the spin sector 
the excited states are separated by the energy uja ~ 
M. So for temperatures T <C uJp,uJa the contribu- 
tion of the excited states to the partition function 
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and, hence, to the specific heat can be neglected. 

Now we can calculate the partition function of the 
system taking into account that impurities are con- 
sidered as independent, and treating the phases <j)i 
at the impurity positions as a randomly distributed 
variable. Then the specific heat of the sample with 
many impurities is related to the partition function 
Zi of a single impurity as 



ZTT 



F, = -T\nZ,, 



where riimp is the impurity density. 

After simple calculations we find the contribution 
of a single impurity to the specific heat 

4e" rS/i^e" + (u + h)^e^ + {u - h)^e-''] 

(l-|-4e"cosh/i)^ 

with u = {Ep Q — Ep i — E^ i)/T and h = iibH/T. 

We are interested in the magnetic field depen- 
dence of the specific heat at low temperatures, T ^ 
M. Then typical values of u for different impurities 
are large, |m| S> 1, and the main contribution to the 
magnetic field dependent part of the specific heat is 
given by values of (/>i for which u > 0. This is natu- 
ral because only those impurities can contribute to 
a magnetic field dependence of the specific heat for 
which the states with a localized spin have the ener- 
gies smaller, than the spinless states. Finally, we find 
for the magnetic field-dependent part of the specific 
heat 



Ch — bn 



cosh n 



(15) 



where b is the fraction of the impurities with u > 0. 
Note that the shape of Cnih) is independent on the 
exact solutions for zero-mode states. In the simple 
model considered here 6 w 1/2, while in general case 
the factor b can be different. 

Magnetic field dependence p3|) qualitatively 
agrees with the experimental data [9110111] . Hys- 
teretic behavior observed in the experiments can 
be understood as well, because one can expect dif- 
ferent metastable distributions of electrons in seg- 
ments between impurities. Application of magnetic 
field and heating/cooling cycles may easily result in 
transitions between such states with close energies. 

In conclusion, we found that strong pinning of 
the LL and partial violation of the spin-charge sep- 
aration by impurities leads to possibility of local- 
ization of electron spins by impurities resulting in 
magnetic field dependence of the specific heat sim- 
ilar to that observed experimentally. We consider 



this agreement as an evidence in favor of formation 
of the LL state in CDW conductors at low tempera- 
tures. Similar magnetic field dependence is possible 
also in the LL with defects without the CDW. 
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